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1. — Lecture One: Tetrahedral Condensed Matter 

Water acts in many ways, from dissolving salt to creating the medium for all life on 
Earth. It is so important and multifaceted that whole books can be and are written about 
it. See, for example, David Eisenberg's and Walter Kauzmann's timeless and recently 
re-issued monograph on the physical properties of water [1]. Our charge in these three 

lectures is to describe something about this material from a molecular perspective. In 
limiting scope, we focus on thermal fluctuations and their consequences on solvation and 
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Fig. 1. - Geometry of a water molecule (left) and a typical hydrogen-bond (right). 



self-assembly. The presentation is like that of a textbook chapter, not a comprehensive 
review. We make use of statistical mechanics at the level it is treated in Ref . [2] . While our 
focus is on water, what we say applies to much of liquid matter, where good background 
is found in Ref. [3]. 

The perspective we adopt is influenced by the results of computer simulations because 
this approach provides unambiguous microscopic detail, albeit for idealized models. Com- 
bined with experiments and theory, simulation is central to all modern understanding 
of condensed matter. Water is a most important example. Computer simulations of 
water were pioneered in the 1970s by Aneesur Rahman and Frank Stillinger. Many sub- 
sequent advances have validated their general approach and enhanced our understanding 
of water. Reviews of that early work in Refs. [4] and [5] remain informative to this day. 

Our first lecture covers properties of pure water, particularly distribution functions 
related to local arrangements of water molecules. The second lecture is about free energies 
of solvation and how these free energies are related to the statistics of spontaneous 
molecular fluctuations in water. The third and last lecture builds from that stage to 
treat forces of self assembly, especially hydrophobic forces, which act on molecular and 
supramolecular scales. 

We begin at the smallest length scales, those of one water molecule. 

I'l. Molecular Structure. - Figure 1 illustrates the water molecule and its most signif- 
icant interaction — the hydrogen bond. Though the molecule is quite polar, its electron 
density is dominated by the electrons of the oxygen atom. As such, the space-filling vol- 
ume of a water molecule is approximately spherical with van der Waals radius of 1.4 A, 
like that of its isoelectronic partner, neon. Because this volume is roughly spherical, it 
is often convenient to identify the position of a water molecule with the position of its 
oxygen nucleus. The OH chemical bond is about 1 A long and the HOH angle is about 
108°. 
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The hydrogen bond interaction is largely electrostatic in origin, and it is strongest 
when the intermolecular separation, R in Fig. 1, is about 3.0 ± 0.2 A, and the angle, 
6*, is < 20°. This linear hydrogen-bond has a maximum adhesive strength of about 
6kcal/mol, which coincides with lOfceT at room temperature (i.e., 25° C). Such a large 
attraction between molecules is unusual. But the strength is only this high over a limited 
range of relative positions and orientations, so that most of the hydrogen-bond strength 
is lost for 9 > 30° and R > 3.5 A, and the interaction decays as 1/R^ to negligible values 
for separations that are a few Angstroms larger. 

Liquids in general are characterized by nearly-balanced competition between energy 
and entropy. When energy dominates, the material is solid, and when entropy dominates, 
the material is gaseous. Liquid water exhibits this balance and is thus ordinary in 
this and many other respects. Its values for density, viscosity and diffusion constant 
are comparable to those of other liquids at ambient conditions (i.e., 25° C and latm 
pressure). But the energy-entropy competition in water is remarkable because its large 
adhesive energy acts over a limited range of configurations. It makes water fragile in the 
sense that its properties have unusually strong dependence with respect to temperature 
and pressure, as discussed below. 

Frozen water is ice, depicted in Fig. 2(a), with the individual molecules arranged 
in nearly perfect tetrahedral order. The resulting hydrogen-bonding network is nearly 
defect-free, with each water hydrogen-bonded to four other waters. Liquid water, on the 
other hand, is bound together by a disordered network of hydrogen-bonds, depicted in 
Fig. 2(b). Hydrogen bonds are constantly breaking and reforming, so that hydrogen- 
bonded partners switch allegiances regularly. 

On average, a water molecule in the liquid will have fewer than the energetically 
optimum four hydrogen-bonding partners. But the number of hydrogen bonds fluctuates, 
and for a tagged water molecule at any one point in time the number typically ranges 
from two to five. On rare occasions, there are even as few as one or none. The precise 
numbers and their probabilities depend upon the definition of a hydrogen bond, and 
conventions for that definition are somewhat arbitrary. For any reasonable definition, 
however, a basin in the potential energy landscape of liquid water coincides with a specific 
hydrogen-bonding network, whereas transition states or saddles separating basins in that 
landscape coincide with breaking hydrogen bonds in the network, and the probability 
distribution for the number of hydrogen bonds that a given molecule experiences is 
unimodal. 

Without theory and simulation, experimental probes of water structure and dynamics 
are often difficult to interpret. This is especially true for spectroscopy. These experiments 
occasionally motivate pictures of water structure different than those described above. 
One example is the two-state model of water, where the liquid is imagined to be an ideal 
mixture of bonded and non-bonded molecules. Unions of experiment with theory and 
simulation have dispensed with that idea, as illustrated, for example, in Refs. [6] , [7] and 
[8]. 

One experimentally-accessible measure of liquid structure is a pair distribution func- 
tion. For water, there are three: 500 (^)) donir) and gimir), which measure the prob- 
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(a) Ice 



(b) Liquid Water 



Fig. 2. - Hydrogen-bonding networks of condensed water, where dashed lines depict hydrogen 
bonds. The picture on the left is a rendering of the molecular structure of ice I. The picture on 
the right is a configuration generated through a classical trajectory (i.e., molecular dynamics 
computer simulations) in which forces are computed according to the SPC/E model of inter- 
molecular potentials [9]. Measures of structure, dynamics and thermodynamics obtained with 
this and similar models are in reasonable harmony with those obtained from experiments on 
real water. 



ability that two different atoms, oxygen or hydrogen, are a distance r apart. More 
precisely. 



(1) 



where p is the number density of water, V is the volume of the system, r^'^'' denotes 
the position of the oxygen atom of molecule z, the angular brackets denote a thermal 
average, and the (5-functions are Dirac's functions, which have unit volume and are non- 
zero only at points where their arguments are zero. Analogous expressions apply for 
gonir) and gimir). In an isotropic system, they depend on only the magnitude of r, 
r = |r|. 

Neutron and x-ray scattering experiments probe different weighted combinations of 
these functions. By adjusting isotopic concentrations and thereby altering weights, each 
distribution function may be estimated. Figure 3(a) shows the three pair distribution 
functions as deduced from neutron scattering measurements. These functions are fully 
consistent with and supportive of a picture in which liquid water is organized with a 
preference for local tetrahedral ordering. 

For comparison. Fig. 3(b) shows the radial distribution function of argon. Argon has 
no strong and directional adhesive forces, and it is therefore called a non-associated liquid. 
Most liquids are non-associated. Structure in such liquids is dominated by packing, and 
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Fig. 3. - Left: Pair distribution functions of water at ambient conditions, as obtained from 
neutron scattering experiments (adapted from Ref. [10]). The shortest distance peaks for the 
OH and HH functions refer to intramolecular separations between atoms. Neutron scattering 
detect these peaks along with intermolecular features. The intramolecular peaks are shown here 
for purposes of comparing with remaining intermolecular features. Right: Radial distribution 
function of a model of Argon at temperature T = 94.4° K and density p = 1.374g/cc (adapted 
from Ref. [11]). The schematic pictures at the bottom depict the pertinent geometries for water 
(left) and for argon (right) from which the main features in the exhibited distribution functions 
can be rationalized. 



since argon atoms are spherical, liquid argon structure is essentially that of a dense fluid 
of hard spheres. An integral over the first peak in g{r) indicates that an atom in liquid 
argon has on average about 12 nearest neighbors, while the corresponding integral over 
9oo{f) indicates that a water molecule in liquid water has on average about 4 nearest 
neighbors. Further, while the second peak for the g{r) of liquid argon is at about twice 
the nearest-neighbor separation, the second peak oi goo{r) is located at around 1.6 times 
the nearest-neighbor separation, corresponding to water molecules that share a common 
hydrogen-bonding partner with local tetrahedral order. 
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Fig. 4. — Pair distribution functions of water at 25° C and 75° C, determined from X-Ray scat- 
tering measurements. Adapted from Ref. [13] (data smoothed for clarity). 



Because non-associated liquids are dominated by packing forces, as detailed in Weeks- 
Chandler- Andersen theory [12], the structures of these fluids are nearly athermal, re- 
sponding to changes in temperature only insofar as these cause changes in density. 
The structure of water, on the other hand, is very sensitive to temperature. Figure 4 
shows the function .9oo('') for water, as determined from X-ray scattering experiments 
at 25° C and 75° C. The most striking feature is that most of the structure beyond the 
first peak disappears above around 50° C, even though the density of the liquid decreases 
by less than 1%. This notable change can be rationalized in terms of the competition 
between entropy and enthalpy. While hydrogen-bonds are energetically very favorable, 
the open structure they impose is entropically unfavorable. At ambient conditions, en- 
thalpy narrowly dominates, while at higher temperatures, entropy plays the larger role. 
In warm water, above 50° C, hydrogen bonds are often broken and no longer dictate lo- 
cal tetrahedral ordering of molecules. This sensitivity to temperature implies that water 
has a large heat capacity at ambient conditions. Indeed, per water molecule, the heat 
capacity at constant volume, Cv, is around 10 fee- For simple non-associated liquids, the 
corresponding heat capacity is an order of magnitude smaller. 

The competition between hydrogen-bond enthalpy and packing entropy is manifest in 
many other properties of water around ambient conditions, though the specific crossover 
temperatures and pressures vary slightly with different observables. For example, at 
atmospheric pressure, water attains a maximum density at 4°C, where its structure is 
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most open (Fig. 5(a)). Other examples are discussed below. 

1'2. Diffusion. - Fluctuations in the structure of liquid water allow molecules to dif- 
fuse. The temperature and pressure dependence of the self-diffusion constant of water 
reflects the degree to which fluctuations change with external conditions. For temper- 
atures below about 310 K, there exists a pressure at which the diffusion constant D of 
water is maximal, as shown in Fig. 5(b). As with the density maximum, this maximum is 
a manifestation of the nearly-balanced competition between enthalpy and entropy, about 
which we will have more to say. First, however, note that the diffusion constant of water 
under most common conditions is, to within an order magnitude, equal to 10~^ cm^/s, a 
value that is typical of dense liquids at ambient conditions. This constant measures the 
rate at which the mean-squared displacement of a particle increases with time, 

(2) {\rf\t) - rl°'(0)|2) ~ 6Dt, for t large. 



The value of D for water is such that after about 1 ps, on average a water oxygen travels 

about 1 A. By taking the timc-dcrivativc of both sides of the above equation, wc obtain 
the relationship connecting the self-diffusion constant to the velocity auto-correlation 
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function, 



(3) 



/3mD = (v^)-^ / dt (v(0) • v(i)) = r„ w 10 



Jo 



Here, /3~^ = UbT, m is the mass of an oxygen atom, and v is shorthand for vj^ , the 
velocity of the oxygen atom of water molecule 1. The right-hand side of this equation 
is the time-scale t„ beyond which (|r[°^(t) — r^°^(0)p) goes as dDt, and this time for 
water, of the order of 1 ps, is also typical of most dense liquids at ambient conditions. 

The diffusion constant is a measure of average molecular motion. The actual motions 
of liquid water molecules are intermittent, and at any single 1 ps time frame, a tagged 
water may or may not move a distance of the order of 1 A. For a water molecule to 
diffuse, it must break an existing hydrogen bond and form a new bond. Conversely, after 
a hydrogen bond is broken, the possibility of the previously-paired water molecules to 
drift apart is much larger than when they are bonded. Diffusion and hydrogen-bond 
dynamics are thus coupled. At least in part because of this coupling, hydrogen-bond 
lifetimes have a non-trivial distribution Phb(t), whose mean corresponds to the typical 
bond lifetime. 



Direct experimental measures of P]m{T) are not available. Its general features can 
be inferred from simulation. The typical lifetime thb is comparable to the diffusion 
timescale Ty. 

The maximum in the diffusion constant reflects the effect of pressure on the entropy- 
enthalpy balance of the hydrogen-bonding network of water. As pressure increases from 
1 bar to about 1 kbar, hydrogen bonds are destabilized, so it is easier for water molecules 
to break and reform hydrogen bonds and thus, diffuse. For moderate pressures, this 
effect is larger than the counteracting slowdown in diffusion expected from the higher 
densities induced by higher pressures. There is a crossover at around 1 kbar, above 
which the diffusion constant, as in a simple liquid, decreases with increasing pressure. 
At high temperatures, as evidenced in Fig. 4, the hydrogen-bonding network is already 
substantially perturbed, and only the simple-liquid dependence of diffusivity on pressure 
is observed. 

1'3. Chemistry in water. - Water plays important roles as the medium or solvent 
for chemical processes. It also plays important roles as a reactant. Chemical proper- 
ties of liquid water in those cases mostly pertain to the presence of disassociated water 
molecules, the hydroxide and hydronium ions, HO" and HaO"*", respectively. They are 
in dynamic equilibrium with intact water molecules: 2H2O ^ HO" -|- H3O"'". While 
the concentrations of these ions are extremely low (about 10~^ times smaller than that 
of intact water molecules), the presence of these ions is crucial to acid-base chemistry, 
biological pumps and motors, and, possibly, transportation fuels of the future. 



(4) 
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A proton in water, sometimes denoted by [H+]aq, is usually bonded to or strongly 
associated with one or more water molecules. Most often it appears in a hydronium ion, 
HaO"*", surrounded by other water molecules. Such protons diffuse in water faster than do 
the individual molecules. They do so by moving along hydrogen-bonding wires, whereby 
a hydronium ion (H3O+ = [H+]aq) that is hydrogen-bonded to a water molecule (H2O) 
transfers the extra proton along a hydrogen bond. Recall that in a hydrogen bond, the 
chemical 0-H bond on one molecule points towards the O of another molecule, OH • • • O. 
The H in that arrangement is called the donor hydrogen, and it is a proton of a donor 
hydrogen that moves from the H3O+ to the adjacent H2O. After this transfer, the 
initially donating ion becomes a stable water molecule and the initially receiving water 
molecule becomes an ion. The new ion can then transfer one of its three protons in a 
similar fashion to yet another water molecule. This chain of events, called the Grotthus 
mechanism, can continue indefinitely provided water molecules are networked together 
with hydrogen bonds. Through this mechanism, excess charge of the hydronium ion can 
travel with minimal diffusion of water molecules. 

While the diffusion of the proton is faster than that of intact water molecules, it is 
only so by factors that arc of order one. This is because the shuttling of charge along 
a chain of hydrogen bonds requires forces from molecules surrounding the chain, and 
these in turn require some changes in the orientations of water molecules — the same 
sorts of reorganization that can also lead to hydrogen bond breaking and molecular dif- 
fusion. Illustrative computer simulations of this dynamics of water and solvated protons 
in Refs. [16], [17] and [18] employ models in which species can break and make chemical 
bonds. 

Auto-ionization of water, [H20]aq — >■ [H+]aq + [OH^]aq, also involves the Grotthus 
mechanism, as shown in Fig. 8 below. But before describing how this ionization occurs, 
we consider further the sorts of fluctuations that occur in water. It is fluctuations — in 
this case, fluctuations of electric fields — that can cause a proton to leave one oxygen and 
join another. 

1'4. Density fluctuations. - Microscopic fluctuations in density are related to how 
water adjusts to the presence of solutes in general, and these fluc;tuations control the 
organization of apolar molecules in water in particular. Remarkably, density fluctuations 
on small length scales obey Gaussian statistics to good accuracy. We will exploit this 
fact in the next two lectures. 

One way to study density fluctuations is to use computer simulations to calculate the 
probability Pv{N) for observing TV waters within a probe volume v. Figure 6(a) shows 
these probabilities for water in spherical probe volumes, and compares them to Gaussian 
distributions with the same mean and variance. Gaussian statistics are not unique to 
water. Figure 6(b) shows a similar collection of Pv{N) distributions in a hard-sphere 
fluid. Gaussian statistics of this sort are consistently found in numerical simulations of 
a wide variety of dense homogeneous fluids. This result is remarkable because it is not 
obvious why the central limit theorem should apply for molecular-scale probe volumes v. 

For large probe volumes v, mean-square density fluctuations are related to the isother- 
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Fig. 6. - (a) Probability of seeing N waters in a splicrical volume v of radius (left to right) 
2.5 A, 3.75 A and 5.0 A. Simulation results at ambient conditions (symbols) are compared with 
a Gaussian of equal mean and variance (lines). After Ref. [19]. (b) Same, for a fluid of hard 
spheres of diameter a at density p = 0.5 cr-^, with a spherical probe volume v of radius (left to 
right) a, 1.5a and 2a. Adapted from Ref. [20]. 



mal compressibility. In particular, 

(5) mn, dp^ 

where p is the average (bulk) molecular density, p is pressure, and 6N = N — {N)y is 
the fluctuation in the number of molecules in v, N. Equation (5) is a standard result of 
the grand canonical ensemble; sec, for instance, Refs. [2] and [21]. The limiting value of 
Eq. 5, dp/d(3p, is a unit-less characterization of the compressibility. It is 1 for an ideal 
gas, and it far exceeds 1 and ultimately diverges as a fluid approaches a critical point. 
But liquid water at standard conditions is not an ideal gas and it is far from a critical 
point. Rather, it is near its triple point. For typical liquids near triple points, dp/d/3p 
is of the order of lO"-^. For water, its value is about 0.06, reflecting the not atypical but 
relatively open and malleable structure of this liquid. 

From Eq. 5, one expects that the variance of Pv{N) will increase as the size of the 
probe volume v = {N)y/p increases. This behavior is indeed observed in the probability 
functions shown in Fig. 6. But for volumes v smaller than about 1 nm^, the mean- 
square size of fluctuations in N is larger than that estimated from Eq. 5. In other words, 
homogeneous water is generally stiffcr on large length scales (more than 1 nm) than on 
small length scales (less than 1 nm). In the next lecture, we have more to say about this 
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Fig. 7. - (a) Schematic of a solute in either state I or state II and the associated energies of 
interaction with the surrounding water. Photoexcitation promotes a transition from I to II at 
time t — (following Ref. [24]). (b) Probability of observing an energy gap £ a time t after 
photoexcitation (symbols: simulation, lines: Gaussian with equal mean and variance). Adapted 
from Ref. [25]. 



fact, which pertains to the correlation length of water. 

1'5. Electric field fluctuations. - Microscopic fluctuations in electric fields play a 
prominent role in charge transfer reactions. In such processes, charge initially local- 
ized on one site transitions to be localized at a different site. These transferring charges 
can be associated with an electron or a proton. When the transition occurs sponta- 
neously, it is the result of thermal fluctuation. In particular, these fluctuations can 
produce electric potentials that make the two localized charge states isoenergetic, facil- 
itating movement of charge from one state to the other. The activation energy for this 
process is therefore dominated by the energies of the facilitating fluctuation in electric 
potentials, which ultimately refers to the free energies and therefore probabilities for 
fluctuations in arrangements of water molecules. These ideas form the basis of Rudolph 
Marcus' celebrated theory of electron transfer. An extended introduction to that theory 
is found in Ref. [22], and in Marcus' Nobel Prize lecture published in Ref. [23]. 

The basic idea of Marcus theory is shown in Fig. 7(a). The potential energy of two 
charge states, denoted I and II, is shown as a function of the structure of the solvent. 
Generically, given a particular solvent configuration, the energies of the two charged states 
will differ by an amount £, called the energy gap. Without an external perturbation, 
the activation barrier to charge transfer is given by the free energy to change the solvent 
structure from a typical configuration to one where the two potential energy surfaces 
intersect. 
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In experiments, rather than wait for charge transfer to happen spontaneously (i.e., in 
the dark), one can inject a photon of energy £{0) to produce a transition at time t = 
from electronic state I to electronic state II. In computer simulations, one can take a 
typical equilibrium configuration of state I and instantaneously change the charge con- 
figuration to coincide with that of state II. In both cases, one can subsequently monitor 
the relaxation of the surrounding solvent as it adapts to the presence of the new charge 
state. Figure 7(b) illustrates the statistics of the energy gap as a function of time elapsed 
since the charge transfer reaction. Here too, Gaussian statistics prevails, in this case for 
the energy gap, not only at equilibrium, but at all times during the non-equilibrium 
relaxation after the charge transfer. 

Another important observation illustrated in Fig. 7(b) is the speed at which the excess 
energy added to the system at t = is dissipated into the solvent. In this example, 
about 60 ksT is injected into the system. Within 30 fs, two thirds of that energy has 
been transferred to the solvent. The fast dynamics illustrated in Fig. 7 are polarization 
fluctuations associated with small rotational or rocking motions of water molecules — 
termed librations. These motions occur with high frequencies because protons have 
small mass and thus water molecules have small moments of inertia. 

Because librations are so rapid, nuclear quantum effects can be important. Any mo- 
tion of frequency ui such that hw/k-oT > 1 will show quantization effects, and at ambient 
conditions, h/ ksT « 25 fs. Indeed, one mechanism of polarization fluctuations in water 
is through quantum tunneling. Librations move molecules through classically forbidden 
regions of configuration space. If the hydrogens are replaced by deuteriums, tunneling 
probabilities decrease because tunneling probabilities in general decrease exponentially 
with increasing particle mass. This decrease explains some of the substantial isotope 
effects in electron transfer reactions. Other observable nuclear quantum effects include 
a small difference in melting temperatures (4°C for D2O, 0°C for H2O) and the tem- 
perature of maximum density (12°C for D2O, 4°C for H2O). Theory for such effects is 
described in Refs. [26] and [27]. 

1'6. Water auto-ionization. - To conclude this first lecture, we discuss the role of 
fluctuations on the mechanism of water auto-ionization. At equilibrium in water, once in 
roughly 10 hours, a tagged intact water molecule will dissociate to become a hydroxide 
(OH^) ion, thereby giving up one of its protons to the surrounding liquid. The process is 
both rare (occurring only once in hours) and fleeting (completing in only picoseconds after 
it starts). Such events can be examined in compiiter simulations through an importance 
sampling of trajectory space. This sampling was used to harvest many independent 
examples of auto-ionization events in water, and one example is shown in Fig. 8. Prom 
analyzing the harvested trajectories and characterizing their transition states, it has 
been found that auto-ionization follows from the coincidence of two fluctuations, which 
we describe now. 

First, a large electric fleld fluctuation briefly destabilizes an 0-H chemical bond and 

thus stabilizes an ion-separated configiiration. Fluctuations large enough to break chem- 
ical bonds occur in water, but they are extremely rare. They are unusual circumstances 
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where electric (mostly dipolar) fields from many, many molecules focus at the site of 
one 0-H bond. But when this happens, because molecular orientations change quickly, 
the large destabilizing electric field lasts for no more than a few tens of femtoseconds. 
The presence of such large and fleeting electric fields in water was discovered by Gra- 
ham Fleming and his co-workers [24] . For the brief period of time when this large field 
persists, an 0-H chemical bond is destabilized and charge separation can occur along a 
hydrogen bond wire by the Grotthus mechanism, as seen in Panels B, C and D of Fig. 8. 

A second fluctuation is required to break the wire of hydrogen bonds along which the 
charge separation proceeds by the Grotthus mechanism. Otherwise, when the rare electric 
field disappears, the separated charges will quickly recombine after the rare electric field 
disappears. The breakage of the wire is evident from Panels D and E of Fig. 8. It 
removes the direct pathway to ion recombination. Average hydrogen-bond lifetimes are 
of the order of picoseconds, while the destabilizing electric field persist for only tens of 
femtoseconds. So, it is unusual for a breaking of a hydrogen bond wire to occur during 
the specific period of time when the electric field acts on that particular wire. After 
the wire is broken, with the charges separated, and rapid recombination is inhibited, the 
solvated proton (i.e., the hydronium ion) can then diffuse away from its parent hydroxide. 

According to these theoretical results, therefore, auto-ionization is the coincidence of 
a large but fleeting electric field fluctuation and a switch in hydrogen bond allegiance. 
Both are rare events, the former much rarer than the latter, and the two must occur 
simultaneously — at the same point in space and time. Figures 8 D and E show typical 
configurations just before and after the system crosses the transition state for this process. 
Experimental observations find that the recombination of initially separated hydroxide 
and hydronium ions is diffusion-limited with a reactive inter-ionic separation of about 
8 A. The mechanism exhibited in Fig. 8 explains why this length is significantly larger 
than the typical separations between nearest-neighbor and second-nearest-neighbor water 
molecules. But direct experimental tests of this mechanism for the fundamental kinetic 
step of pH are not yet available. 

2. — Lecture Two: Solvation 

The previous lecture describes the behavior of pure liquid water, stressing how its 
local structure is tetrahedral, and how its density and polarization fluctuations obey 
Gaussian statistics to a good approximation. This second lecture describes how the na- 
ture of these fluctuations determines water's behavior as a solvent. The central quantity 
to be considered is solvation free energy — the reversible work done on the solvent to 
accommodate a solute molecule. This quantity determines the probability of solvation 
and its associated driving forces. 

These forces can be strong, such as when water successfully outcompetes powerful 
ionic bonds, or when it completely shuns solutes, as happens in oil-water de-mixing. 
Amphiphilic substances add geometric frustration to the physics of solvation, where each 
molecule has two distinct portions at a fixed separation, one that binds water and the 
other that repels water. The fixed separation is a geometric constraint that frustrates 




Fig. 8. - Typical 150 fs trajectory illustrating an auto-ionization event, from Ref. [28]. Each 
panel is separated by 30 fs. The transition state for the process is visited in a time frame 
between those pictured in Panels D and E. The hydroxide and hydronium ions are shown in 
blue and yellow, respectively. The subset of hydrogen bonds that play a specific role in this 
particular trajectory are shown with dotted lines. 



water structure, which in turn can lead to aggregated nano structures, hke membranes 
and micelles. Proteins are large amphiphilic molecules for which water significantly 
affects dynamics and thermodynamics, for instance by providing one of the driving forces 
for protein folding and assembly. 

This lecture describes underlying principles and simple quantitative estimates of free 
energies related to these behaviors. A few consequences are described in the next lecture. 

21. Solvation free energies. - Imagine a system where initially there is only bulk 
water and a solute in vacuum far away. The solute is then slowly inserted into the water 
(Fig. 9). This insertion requires performing some reversible work, which is equal to the 
free energy difference A/i between the system before and after insertion. This free energy 
can be calculated from statistical mechanics as a ratio of partition functions, 

W Jdxe-f'^oi^) ' 



where x denotes a configuration of the system, Eq{x) is the energy of the pure solvent and 
Ei{x) is the energy of the solvent when the solute is in the liquid. The energy function 
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water water + solute 1 water + solute 2 



Fig. 9. - Solvation free energies and the forces of assembly. The excess chemical potentials for 
the separated pair (the red and green particles in the middle picture) and the complexed pair 
(the red and green particles in the right-hand picture) are the solvation free energies A/ii and 
A/i2, respectively. The change in solvation free energy, A/i2 — A^i is the reversible work done 
on the liquid to assemble the complexed red-green solute from the separated red-green pair. 



Ei{x) is parameterized by the state of the solute, so changing the condition of the solute 
generally changes that function. 

The quantity A/i is called the solvation free energy. In a dilute solution, where the 
solutes do not interact with each other, the chemical potential /x of the solute is given by 

(7) /3M = /3A/i + ln(p/po), 

where p is the concentration of the solute, relative to its concentration po in the standard 
state. The logarithm on the right-hand-side is the only contribution to when the solute 
and the water do not interact. This logarithm is thus the ideal gas chemical potential. 
For this reason, the solvation free energy is sometimes also called the solute's excess 
chemical potential — the excess relative to that of an ideal gas. One way to deduce excess 
chemical potentials from experiments is to measure partitioning fractions of solutes in 
two coexisting fluids. Connections like that between A/x and experiments are discussed 
in standard texts, e.g. [2] and [29]. 

When two or more solutes are dissolved in water, their solvation free energy can 
depend on their relative positions. This dependence may be due to packing effects, 
for example, or the structure of the hydrogen-bonding network of water around one 
solute being perturbed by the presence of the other solute. Regardless of how this 
dependence arises, the solvation free energy of the two solutes when they are nearby 
(A/i2 in Fig. 9) may differ from its value when they are far apart (A/j,i in Fig. 9). If the 
difference A/X2 — A/^i is negative, then the associated state is favored. If it is positive, 
on the other hand, the dissociated state is favored. The solvent provides a force of 
assembly or disassembly, respectively, depending upon whether A/i2 — A/ii is negative or 
positive. This type of force is what drives oily particles together in water and is ultimately 
responsible for the hydrophobic effect, the main topic we discuss in Lecture Three. 

The excess chemical potential is but one particular free energy that can be computed 
for a solute-water system. In particular, consider a variable A that interpolates between 
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one solute and another. For instance, A may be the distance from the solute to the liquid- 
vapor interface, or it may be the distance between two dissolved solutes, or a parameter 
that is used to progressively "create" or "grow" solutes into the solvent. Denote by 
A(x) the value of A for a particular configuration x. At equilibrium, configurations with 
specific values of A occur with probability 

(8) P(A)(x j dxSiX{x)- X)exp[-PEi{x;Xix))]. 

The right-hand side of this equation is a partition function for configurations constrained 

to the surface A(.t) = A. The quantity — fcBTlnP(A) is therefore the free energy of the 
system constrained to have that particular value of A. When A = 1 and A = correspond 
to the solvent with and without the solute, respectively, then —kBTln[P{l)/P{0)] = Afj,. 
At intermediate values of A, — fcBTln[P(A)/P(0)] serves as a useful generalization of 
solvation free energy. 

We now apply this generalization in two ways, first to estimate solvation free energy of 
excluded volume (i.e., the reversible work to make space for a solute in a liquid solvent), 
and then to estimate solvation free energy due to the charge of a solute (i.e., the reversible 
work to polarize the solvent). 

2'2. Solvation of small excluded volumes. - To compute the solvation free energy for 
an excluded volume v, we use for A(a;) the observable N{x), which denotes the number 
of water molecules in a volume of that size in the liquid. In other words, we consider the 
probability that A'' molecules exist in a volume v, Py{N). The excess chemical potential 
for a particle of volume v absent any other interactions with the solvent is the reversible 
work to empty that volume, which is to say 

(9) A;U^ = -fcBTlnP„(0). 

Gerhard Hummer, Lawrence Pratt and their co-workers introduced this way of consid- 
ering solvation [19], and we follow their lead. 

Specifically, from the empirical observation discussed in the previous lecture — that 
density fluctuations in small volumes are Gaussian — we write 



(10) P,(A)«-=^exp 

V Zn Gy 



2ay 



where {N)y = pv is the average number of water molecules in the volume v, p being the 
mean bulk density of the liquid, and 



(11) 



= {{SNn 
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is the mean-square fluctuation in the number of molecules in that volume. Accordingly, 
the solvation free energy of an excluded volume v is 

(12) A/x„ « kBT 

Equation (12) is an important result that holds to the extent that the Gaussian 
distribution, Eq. (10), is valid. This in turn requires that v not be too large. In particular, 
water at standard conditions is close to coexistence with its vapor phase, and as such, 
emptying a large enough volume can nucleate a vapor bubble and its associated liquid- 
vapor interface. The Gaussian approximation to Pv{N) does not include this physics of 
phase change, but rather includes only the physics of fluctuations that do not move the 
region within v far from its tj^pical liquid state. We will see that this criterion limits the 
usefulness of Eq. (12) to w < 1 nm"^. 

Keeping to v < 1 nm^, Eq. (12) provides a means for estimating solvation free energies 
associated with excluded volume in water. The values of {N)y and ay can be estimated 
simply. Specifically, the average number of water molecules in v is {N)y = pv, and the 
mean-square fluctuations are given by 

(13) ay= [ dr / dr' (5p(r)Jp(r')) , 

where Sp{r) is the fluctuation (i.e., deviation from the mean) of the molecular density at 
position r. The integrand is related to the pair distribution function 500(1"), 

(14) {dp{r) 5p{r')) - pS{r - r') + p'[9oo{\r - r'|) - 1], 

and can thus be calculated when goo{r) is known. Solvation free energies estimated 
from these formulas agree well with those obtained from computer simulations, and 
these formulas form a convenient basis for accurate estimates of hydration free energies 
for small oily molecules in water. The extent of accuracy is illustrated later in this lecture 
with Fig. 14 below. 

Equation (12) shows that the free energy to exclude small volumes is primarily en- 
tropic. In particular, this free energy grows with increasing temperature. To the extent 
that excluded volume effects are dominant, therefore, small solutes become less soluble, 
not more, as temperature increases. This trend is indeed observed in experimentally 
measured solubilities of small apolar molecules in water. A trove of such data is found 
in Charles Tanford's classic monograph [30]. At temperatures above about 50° C, how- 
ever, the structure of water begins to change substantially (Fig. 4), so that the solvation 
free energy of a small solute stops being proportional to temperature. The theoretical 
predictions are graphed in Fig. 10. 

To gain hirthcr insight, consider v significantly larger than a correlation volume. Such 
sizes are at the upper limit of where the Gaussian approximation is accurate. But for 
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Fig. 10. - Predictions from Eq. (12) of solvation free energies for hard spheres in water, with 
hard-sphere radii comparable to those of Ne, Ar, Me (methane) and Xe, as. Around T w 400 K, 
all curves have the same slope (i.e., the same solvation entropy). Adapted from Ref. [31]. 



such volumes, the compressibiUty theorem mentioned in the previous lecture gives 



(15) 



(pv) ^ a-u 



{{SNf), dp _ 



An approximate expression for the solvation free energy is thus 



(16) 



A/i w kBT 



pv 
2X 



1 



\n{2npvx) 



This simple estimate is quahtativcly correct, but overestimates the more generally ac- 
curate Eq. (12). Nevertheless, it immediately shows that hydration free energies due to 
excluded volume of small apolar species are typically 5 to 10 ksT because the quantity 
X has a value of about 0.06 for water for a range of temperatures around ambient con- 
ditions. It further illustrates that the solvation free energy grows roughly linearly with 
increasing volume in the regime where the Gaussian approximation is valid. 

Figure 11 shows how the experimentally measured solvation free energies of small-to- 
medium alkyl chains scale with chain length. These solutes are extended but thin, so 
they can be treated with the theory we described above for small solutes. As chain length 
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Fig. 11. - Solvation free energies of CH3(CH2)n-iCOOH as a function of chain length. The line 
is a linear fit to the data. Adapted from Ref . [35] . 



and solvent-excluded volume are proportional, these results iUustrate that solvation free 
energies of small solutes are proportional to excluded volume. In many discussions, 
e.g., Ref. [32], these same results have been used to argue that solvation free energies 
of even small solutes may be modeled as scaling with solvent-accessible surface area 
(SASA), because chain length and SASA are also proportional. However, the fact that the 
apparent surface tension that would result is much smaller than the measured water-oil 
surface tension, and that the alkyl chains become less soluble, not more, as temperature 
is increased, shows that the proportionality with SASA is a misleading coincidence. 
References [33] and [34] provide further discussion on this point. 

Equation (12) rationalizes a phenomenon called entropy convergence. To discuss it, 
we introduce the entropy of solvation, which is the entropic component of the solvation 
free energy, given by 

It is found that along the saturation curve of water, the entropy of solvation of many 
different small solutes is different, but surprisingly, converges to a common, small value 
at a temperature of about 400 K. This convergence is illustrated in Fig. 10. Around 

this temperature, it is found that the factor is nearly athermal, so the solvation free 
energies are essentially proportional to Tp^ (T) . This combination is non-monotonic with 
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(a) (b) 

Fig. 12. - (a) Calculating the solvation free energy of an ion in two steps. A free energy A/i„ 
is needed to expel waters from the ion's hard core. A free energy A^, is then required to 
place a charge q in the ion and thus polarize the solvent, (b) Free energies of the potential <j) 
at the center of an ion, before (solid) and after (dashes) the charge is inserted. The induced 
polarization of the solvent shifts the mean potential to = —q/k.. 

temperature, and is maximal around T — 400 K. Thus, its derivative is zero there, leading 
to a vanishing entropy of solvation at that temperature. The observed small but non-zero 
value arises from the second term of Eq. (12) and the small temperature dependence of 

A similar convergence is also found around T — 400 K in the per-residue entropies of 
unfolding of proteins. Though it superficially resembles the entropy convergence of small 
solutes, the origin of this effect is actually more complicated and due to the interplay of 
small and large length-scale physics that we discuss in general at the end of this lecture 
and in Lecture Three. Reference [36] provides further discussion on this particular point. 

2'3. Solvation of ions. - So far, we have discussed solvation of excluded volume. While 
all solutes exclude water from some region of space, most solutes act with additional forces 
on water, including weak dispersive attractions and electrostatic interactions. Weak 
interactions change the solvation free energy quantitatively but not qualitatively. On 
the other hand, owing to the large dielectric constant of water (e « 80), electrostatic 
interactions are significant, and often dominate solvation behavior when they are present. 
For instance, since we know that water dissolves many salts, it must be the case that the 
solvation free energies of the separate ions is larger than the strength of the ionic bonds 
between them in a crystal, which is about 1 eV, or about 40 fceT at ambient conditions. 

As a model for electrostatic interactions, we estimate the solvation free energy of a 
single ion, as outlined in Fig. 12(a). An ion typically interacts with water by excluding a 
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water from its hard core and by electrostatically polarizing the surrounding solvent. To 
be concrete, a Cl~ ion excludes water from a sphere with radius of about 2.5 A, and can 
be modeled as having a point charge of — e at its center. The free energy A/z^, to expel 
the water from the hard core is given by Eq. (12). The free energy A/Xg to then place a 
point charge q at the ion is given by 

rd(ie-''[^°('^)+«"^l 
^ ' /d<^e-/3^o(« ' 

where 4> is the value of the electrostatic potential at the center of the ion given that a 
cavity for it has been created, and FQ{(f>) is the free energy associated with a particular 
value of this potential, equal to —kBT\nP{(j)). As outhned previously, the statistics of 
electric field and potentials in water are Gaussian, so Fo{(p) is quadratic in (p. Moreover, 
since the water is charge neutral, the mean value of ^ should be essentially zero(^). 
Hence, we approximate 

(19) Fo(<^) « ^ncp\ 

where k is related to the fluctuations of (j) in the absence of the charge by equipartition, 

(20) mr) = ^. 

These fluctuations are, in turn, controlled by the dielectric constant of water and the 
size of the ion, as discussed below. With this expression for Fo{(l)), we can evaluate the 
right-hand side of Eq. (18) and simplify the result to obtain 

(21) A„, = -l(3q'mf} = -^. 

When the charge is placed into the ion, the surrounding solvent is polarized. As 
depicted in Fig. 12(b), the potential fluctuations are unchanged but the mean value of 
the potential shifts to 

(22) (</.), = -q/K 



(^) It wont be exactly zero, however, in the presence of large enough solutes that cause signifi- 
cant inhomogeneity. Liquid-vapor interfaces have a small surface dipole moment density, since 
a small number of water molecules at the interface tend to align their dipoles with the surface 
normal. The bulk liquid has been measured to be at a potential of about 0.1 V higher than 
the bulk vapor. This potential difference is the so-called zeta potential of the interface. For a 
discussion of these effects and how the presence of ions modifies the structure of the liquid- vapor 
interface, see Ref. [37]. 
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This induced potential acts as a reaction field of the solvent to the presence of the charge. 
Polarizing the solvent in this way has a free energy cost of 

(23) Fom,) = I'^iQ/'^f = 1^, 
but results in a favorable interaction energy of 

(24) = 

between the ion and the polarized solvent. 

The value of k can be determined by relating the fluctuations in the dipole moment 
density field of water to its dielectric constant and the geometry of the cavity (for details, 
see Ref. [38]). It is found to be 

(25) k = r(i-^^ \ 

where R is the radius of the ion. This result is identical to the continuum electrostatics 

result [39], where the work required to move a charge q from vacuum into a spherical 
cavity of radius R carved into a dielectric material with dielectric constant e is 

(26) ^^^ = -U'-l)- 

Equation (26) is the Born solvation formula. It is negative, so it is always favorable to 
solvate an ion. It scales as g^, so divalent ions are much more strongly solvated than 
monovalent ones. Finally, it scales as 1/R, so smaller ions are more soluble than larger 
ones. 

For water, e is so high that the term in brackets is virtually unity. To within an order of 
magnitude, the solvation free energies of ions are about e^/(l A) « lOeV ~ 400 /cbT. As 
anticipated, these free energies are indeed sufficiently large to overcome the ionic bonds in 
salt crystals. Figure 13 shows a more detailed comparison of the Born solvation formula 
to experimental enthalpies of solvation. It is clear that the Born solvation formula works 
quite well for all of these cases, and thus captures the dominant physics of solvation for 
ions. 

2'4. Solvation of large solutes. - In both of the scenarios considered above, we could 
calculate solvation free energies from the statistics of small structural fluctuations in 
water. The solute, whether ideal or charged, acted as a small perturbation that did 
not fundamentally change the nature of the liquid. As solutes increase in size beyond 
about 1 nm, however, water responds to their presence in a massively collective fashion, 
undergoing nano-scale reorganization manifesting a macroscopic phase transition — the 
liquid-vapor phase transition — that occurs at thermodynamic conditions very close to 
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Fig. 13. - Enthalpies of hydration (A// — T[9A/x/9T]) for a variety of monovalent ions as a 
function of inverse ionic radius, as measured experimentally (circles) and predicted by Eq. (26) 
(line). Adapted from Ref. [40]. Similar correlations are seen for polyvalent ions. 



those of ambient water. The solvation free energy is then dominated by the surface 
energy of the interface between the two phases, which for macroscopic solutes is given by 

(27) Afi -J> 'yA, (large solutes), 

Here, 7 is the surface tension between water and the solute (about 72mJ/m^ for a solute 
that only excludes volume) and A is the surface area of the solute. 

For nanometer-sized solutes, there is a gradual crossover between the small length- 
scale behavior described by Eq. (12) and the large length-scale behavior described by 
Eq. (27). This crossover is depicted schematically in Fig. 14 for hard sphere solutes of 
volume V. The two estimates for Afi intersect around v of Inm*^. This sets the scale 
for where the crossover is relevant, though corrections to Eq.(12) are already important 
for V > 0.5 nm^. Solvation free energies eventually do tend to jA, but convergence to 
this limit is slow because the ratio of these two quantities tends to 1 only as v~^^^. 
Corrections to Eq. (27) are significant even for excluded volumes measuring tens of cubic 
nanometers. 

The breakdown of Eq. (12) for larger v indicates that the Pv{N) distributions for these 
larger volumes deviate significantly from Gaussian behavior. This is indeed the case. 
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Fig. 14. - Solvation free energies for hard spheres of volume v, computed through direct sim- 
ulation with the SPC/E model of water at ambient conditions (black, Ref. [41]). The small 
length-scale estimate (Eq. (12)) and the large length-scale estimate (Eq. (27)) for the same 
model are also shown. 



Figure 15 shows an explicit calculation of Pt,(iV) where t; is a cube of volume 1.7nm^. 

The breakdown of Gaussian statistics appears as a fat tail in this distribution. This fat 
tail indicates that there is a mechanism beyond Gaussian density fluctuations that can 
be involved when removing N waters from the probe volume v. The fat tail in lnP„(iV) 
is well-described by a function that scales as —^{SN)'^/^, which describes the free energy 
to form an empty cavity of volume SN/p inside of the probe volume. The physical 
consequences of this fat tail and corresponding crossover in length-scale are discussed in 
the next lecture. 

3. — Lecture Three: Hydrophobicity and Self-Assembly 

This lecture is about hydrophobicity, which essentially means it is about how and 
why oily species termed hydrophobic species — separate from water in nanometer scale 
structures. "Hydrophobic" is actually a misnomer, because the underlying physics is not 
about oil fearing water, but rather about hydrogen bonds that are lost when water mixes 
with oil. A hydrophobic species is simply a molecule or complex of molecules that binds 
more weakly to water than water binds to itself. A hard sphere is the simplest (and 
idealized) example. Surfaces of nominally hydrophilic molecules can also be hydropho- 
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Fig. 15. - Water number distribution, Pv{N), for a cube of volume v — 1.7nm^ computed from 
importance sampling with a computer simulation of the SPC/E model of water at ambient 
conditions (solid), and the distribution that would be expected if the statistics of these number 
fluctuations were Gaussian (dashed). The fat tail is an indication of the crossover in solvation 
free energies to the large-length scale regime. Adapted from Ref. [42]. 



bic when the geometry of possible water-surface binding sites is incommensurate with 
favorable hydrogen patterns of liquid water. 

Hydrophobic forces are the water-mediated interactions between oily species in water 
that cause these species to segregate or demix from water. Macroscopic demixing of oil 
and water manifests a first-order phase transition. Such manifestations require clusters 
of size larger than a critical nucleus. A small enough hard sphere in water will not trigger 
the physics of phase separation. Hydrophobic forces therefore become significant only 
for sufficiently large hydrophobic species. We will see that "sufficiently large" implies 
hydrophobic surfaces of low curvature that extend over lengths of the order of 1 nm or 
more. 

Lecture Two provides the background to these ideas. Indeed, the 1 nm scale mentioned 
here is related to a finding discussed in the previous lecture, where we describe how 
solvation free energies are related to statistics of solvent density fluctuations and solvent 
polarization fluctuations. We show there that these statistics are essentially Gaussian for 
cases where the solute does not largely perturb the liquid. But we also show there that 
this situation changes in the disruptive presence of sufficiently extended hydrophobic 
surfaces. This lecture focuses on the consequences of these changes and a theory that 
embodies the underlying physics. Ref. [33] is review of the topic written with the same 
perspective as that of this lecture. 
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Fig. 16. - Driving force AG of hydrophobic assembly. At room temperature (thick black Hues), 
the crossover length scale is about Inm. At higher temperatures (thin red lines), the crossover 
length scale decreases and the driving force increases. After Ref. [33]. 



3'1. The driving force for hydrophobic assembly. - Figure 16 helps to show how the 
assembly of hydrophobic aggregates in water is explained in terms of the length-scale 
dependence of solvation free energies discussed at the end of Lecture Two. As schematized 
in that figure, the solvation free energy to solvate N small, dissociated solutes scales as N, 
and is mostly entropic in origin. The hydrogen-bonding network of water is stretched 
but not disrupted by these solutes. Except for small packing effects common to all 
dense liquids, the total solvation free energy of this collection of solutes is independent 
of their positions. As such, entropy disfavors aggregation, so the solutes disperse evenly 
throughout the solvent. 

If enough of these solutes associate spontaneously, however, it becomes impossible for 
water to wrap its hydrogen-bonding network around the cluster. It then becomes prefer- 
able to sacrifice some hydrogen bonds completely and nucleate a vapor-liquid interface 
around the cluster. The cost of maintaining this interface is primarily enthalpic, and 
scales as N'^/^ . Hence, for large enough N , it is energetically favorable for small solutes 
to associate into a larger cluster, with the free energy difference AG between the dis- 
sociated and associated states termed the "driving force" for hydrophobic assembly. At 
ambient conditions, the critical nucleus that needs to form before aggregation becomes 
favorable is typically about 1 nm in size. 

The temperature dependence of the driving force for assembly is non-trivial. As 
evidenced by Eq. (12), the solvation free energy of very small solutes is proportional 
to temperature, because of its entropic origin. Conversely, the surface tension of most 



Fluctuations in water 



27 




(a) (b) 

Fig. 17. - Micelle assembly, (a) Small amphiphiles dispersed uniformly as monomers in solution 
(three of which are shown); these monomers are in equilibrium with n-mers - clusters of n 
amphiphiles (one of which is shown), (b) For amphiphile concentrations p below pcmc, only 
the monomeric species (solid) is present. Above pcmc, n-mers form (dashes), and adding more 
amphiphiles only increases the number of n-mers 



liquids, including water, decreases with temperature. In other words, whereas rises in 
temperature make small solutes less soluble, they make large solutes more soluble. This 
curious dichotomy in the temperature dependence of solvation free energy has two con- 
sequences. First, the size of the critical nucleus for aggregation is reduced at higher 
temperatures. Second, the driving force towards hydrophobic association is stronger in 
hot water than in cold water. These consequences are pertinent to protein thermody- 
namics, where hydrophobic interactions stabilize the core of most globular proteins. At 
colder temperatures, the driving force for this hydrophobic collapse decreases, so at low 
enough temperatures, proteins undergo cold denaturation. 

At high pressures, the dominant component to solvating large solutes is not the cost 
of forming a liquid-vapor interface, but the work needed to create a vacuum bubble in a 
high-pressure environment, which scales as N . Hence, above pressures of about 500 atm, 
the driving force for hydrophobic assembly disappears, and proteins undergo pressure 
denaturation. 

3'2. Micelle Assembly. - When small amphiphilic molecules, such as the fatty acids 
that compose cell membranes, are dissolved in water at low concentrations, entropy favors 
dispersing them uniformly as monomers. At higher concentrations, however, the driving 
force for hydrophobic assembly overcomes this entropic effect, so collections of these 
molecules assemble into nano-scale clusters, which are called "micelles." In this way, the 
hydrophobic tails of the amphiphiles are separated from water by a layer of hydrophilic 
head groups. Experimentally, the change in aggregation behavior occurs abruptly at a 
critical micelle concentration, yOcmo whose non-trivial temperature dependence can be 
understood from the considerations we have been describing — the competition between 
small length-scale solvation and interface formation. 

A simplified model of micelle assembly is guided by Fig. 17. For simplicity, we shall 



28 



D. Chandler and P. Varilly 




Fig. 18. - Thermodynamic cycle of micelle formation (adapted from Ref. [43]). The free energy 
of micelle assembly AG consists of three contributions: a term AGi for forming a cavity of the 
right size in the solvent, a term AG2 for detaching the head and tail groups and transferring 
the hydrophobic tails from solution into the cavity, and a term AG3 for reattaching the head 
groups to the tails. 

assume that the micelles that form are all the same size, each containing n amphiphiles, 
with that number n to be determined. (For large enough n we expect the dispersion 
of n-mer size, An, to be small as a consequence of a law of large numbers, specifically 
An ^ ^n.) Accordingly, the total amphiphile concentration is 

(28) p = pi+npn, 

where pi is the density of monomers and p„ is the density of n-mers. The law of mass 
action relates these two quantities by 

(29) p„ =prexp(-/3AG), 

where AG is the driving force for assembly of the micelle, namely the free energy differ- 
ence between one n-mer and n monomers. 

To estimate the driving force for assembling an n-mer, we use the thermodynamic 
cycle depicted in Fig. 18. The width of the amphiphile is denoted by a and its typical 
length is denoted by S. The n-mers have radius R and surface area A. The assembly is 
achieved in three steps: 

1. A cavity of the size of a micelle is formed in the solvent. Since this cavity is large, 
the free energy to create it is proportional to its area A, so 

(30) AGi « 7wator-vapor^- 
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In terms of n, a and S, the area A scales as a'^{5/a)'^^^n^^'^ . 

2. The head and tail groups arc detached, at energy cost ebond per molecule, and the 
tails are transferred from the solution to what will be an oily region. This is the 
free energy of transfer A/i^transfor of the tails from water to oil, and depends on 
the tail's size and chemical character. Finally, what was a water-vapor interface 
after Step 1 is now a water-oil interface. Due to dispersive interactions, its surface 
tension is thus reduced by an amount A7 given by 7water- vapor — 7oii- water- Hence, 

(31) AG2 « nCbond - nA/itransfer " A7A. 

3. Finally, the head groups are reattached to the tails, for an energy gain of ebond per 
molecule. The additional entropic cost of restricting the n head groups to be at the 
surface of the micelle is quite large. Since this restriction is exactly analogous to 
the requirement of charge neutrality in a polarized dielectric, it can be estimated 
accurately by analogy [44, 45]. Specifically, the entropic cost of maintaining a 
separation of a head group from a tail group at a distance r goes as 1/r , which makes 
the effects of this cost isomorphic to that of an electrostatic cost. Accordingly, the 
total entropic cost has the same /R functional dependence as the Born solvation 
energy, Eq. (26), where n plays the role of charge and R ~ n^/^ is the micelle 
radius. The pre-factor is related to the aspect ratio of the amphiphile. We thus 
obtain the estimate 

(32) AGs « -nebond + knT{a/5f'''n'''''. 

Above, we focused on the functional dependences of the free energies on n, a and 5. 
Adding up the contributions, we obtain the total free energy of assembly: 

(33) AG « 7oii-watera2(5/a)2/^n2/3 - nA/itrar.sfer + kjiT{a/ Sf/^vP'^ 

For a given set of parameters in the above equation, we find the number of monomers 

n = n* that minimizes AG/n to obtain the typical micelle size, and pcmc can then be 
determined through the law of mass action. We omit the algebra and simply give the 
final result [43]: 

(34) Inpcmca^ ~ C (^7oil-watera^)^^^ - /SA/itransfer, 

with c = (5832/49)1/3 ^ ^ ^ 

gure 19 shows how this prediction compares to the ex- 
perimental measurements. The data can clearly be explained quantitatively by invoking 
only the length-scale dependence of the hydrophobic effect and the geometric constraint 
that head and tail groups be adjacent. 



30 



D. Chandler and P. Varilly 




1.8 



0.8 



8 10 12 14 16 



10 



30 



50 



70 



m 



T (°C) 



(a) 



(b) 



Fig. 19. - Comparison of Eq. (34) with experimental pcmc of CH3(CH2)(m-i) (OCH2CH2)60H. 
The single adjustable parameter a w 3 A is fit to the m = 12 result at T = 25° C. (a) Dependence 
of Pcmc on amphiphilo length at T = 25°C. (b) Dependence of pcmc on temperature for m = 
12. In both plots, the data shown is experimental (circles and triangles, from two different 
experiments) and the result of Eq. (34) (solid). Adapted from Ref. [43] 



3'3. Dewetting transitions in hydrophobic assembly. - As discussed in Section 2'4, and 
as illustrated in Figs. 14 and 15, the length-scale dependence of solvation free energies is 
intimately connected to the formation of liquid- vapor- like interfaces around large solutes. 

Liquid-vapor interfaces are soft modes, which is to say that it costs little free energy 
to deform these structures. Creating a cavity next to a large solute is thus much easier 
than creating the same cavity in bulk water. A demonstration of this effect in water 
is shown in Fig. 20. Here, a model 24 x 24 x 6A plate that is either hydrophobic or 
hydrophilic is placed next to a 24 x 24 x 3 A probe volume v, and the probability Py{N) 
of finding waters in v is calculated. The hydrophobic plate induces a fat tail in the 
PviN) distribution, similar to the one depicted in Fig. 15. This occurs because the 
liquid-vapor interface can be easily deformed to create an empty space inside the probe 
volume. Next to the hydrophilic plate, where the interface between water and the plate 
is not liquid- vapor- like, the Pv{N) distribution is identical to that obtained when the 
probe volume v is far away from the plate. 

The fat tails in Pv{N) distributions are present for large probe volumes in bulk and 
enhanced by a nearby large solute. This behavior has two main physical consequences. 
First, a hydrophobic object that excludes water from a volume v has a lower solvation 
free energy when next to the large solute than when dispersed in bulk water. Indeed, the 
difference in — fcBnnP^,(0) next to the large solute and in bulk measures the reversible 
work needed to move the object from a specific position in the bulk to the vicinity of the 
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Fig. 20. - Effect of a large solute on the density fluctuations of the surrounding solvent, (a) 
The model plate and the probe volume v. (b) Probability P^(N) of observing waters in the 
probe volume v in bulk (solid black), next to the hydrophobic plate (solid blue) and next to 
the hydrophilic plate (dashed red, nearly indistinguishable from solid black). Adapted from 
Ref. [42]. 



large solute. This is precisely the free energy of hydrophobic adhesion. 

Second, these fat tails also signal an underlying phase instability that can be exposed 
by an external perturbation. The idea is illustrated in Fig. 21. Consider an external 
perturbation of the form eN. This perturbation might be the attractive dispersive in- 
teractions between a large solute and the water in its vicinity, or the effective repulsion 
of water from a cavity in confining geometry. If the Py(N) distribution is Gaussian, 
then a linear perturbation results in another Gaussian free energy of equal width but 
different mean number of waters in that volume, (7V)^. Small changes in the strength of 
the perturbation induce small changes in {N)y. If, on the other hand, Pv{N) has a fat 
tail, then a linear repulsive perturbation can result in a precipitous reduction in {N)^. 
This phenomenon is called a dewetting transition. 

This transition can occur in volumes confined by hydrophobic surfaces, such as water 
inside nano-scale tubes. Emptying and filling such volumes is thus collective, with several 
water molecules leaving or entering in bursts. Refs. [46] and [47] illustrate this behavior. 
The phenomenon seems likely relevant to functioning of biological pores. It is akin to a 
liquid-vapor transition, but on a scale of nanometers. Dewetting also appears in cases 
where two hydrophobic protein surfaces approach one another. Water confined by these 
surfaces is destabilized, and as water departs the surfaces move closer together to fill the 
vacated volume. This behavior seems relevant to the dynamics of protein folding and 
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Fig. 21. - Effect of an external perturbation of the form eA'^ on different Pv{N) distributions, 
(a) When v is small and Pv{N) is Gaussian. The perturbation simply shifts the mean number 
of waters (A'^)i, in v. (b) When v is large, or v is next to a hydrophobic solute, Py{N) has fat 
tails. A large enough perturbation results in a first-order microscopic phase transition. In both 
figures, the Pv{N) distribution shown correspond to e = (black/solid), e < (blue/dot-dash) 
and e > (red/dashed). 



assembly, as our group first discussed in Ref. [48]. 

3'4. Theory of Dewetting. - Theory for dewetting requires treatment of both inter- 
faees and small length-scale fluctuations. High free-energetic costs of solvating excluded 
volumes at small length scales gives way to lower costs in the presence of soft interfacial 
fluctuations. Here, we describe how to build a theory that captures this physics with a 
density field that describes interfaces and a coupling of that field to small-length scale 
fluctuations. The development uses some elements of statistical field theory, and while 
it is therefore a step beyond the simplicity adopted in the earlier parts of our lectures, 
good textbooks on the topic do exist. See, for instance, Mehran Kardar's [49]. 

Interfaces are well described by density fields that vary slowly on molecular scales. 
In the simplest case, the energetics of such a field, n(r), is given by a Landau-Ginzburg 
hamiltonian of the form 



(35) 



l3Hi,[n{r)] ^ dr 



where we use subscript "L" to indicate that this hamiltonian applies to a fluid on large 
length scales only. The quantity w{n) is a local (grand canonical) free energy density in 



Fluctuations in water 



33 



units of fee T = 1/^, and the parameter m determines the free energy cost to create an 
inhomogeneity. 

In mean field theory, the average n(r) is the function that minimizes this hamiltonian 

(subject to whatever constraints specify the ensemble considered). This minimization 
produces a spatially invariant value for (n(r)), except at conditions of phase coexistence 
where an interface separates volumes with different values of the field. The interfacial 
tension for this interface is proportional to m and the shape of the interface is determined 
by the function w{n). These relationships can be read about in standard texts, e.g. 



The actual molecular density, p{r), can be written as a slowly varying field like n(r) 
plus a correction, where the correction accounts for fluctuations that occur on small- 
length scales. In particular, we take 



where 6p{r) is the small- length scale part. There is flexibility in defining this decom- 
position, but it is important that n(r) varies little over a length ^, which is correlation 
length of the homogeneous liquid. With this generic criterion, the vapor phase of water 
is where n(r) is close to zero, and the liquid phase is where n(r) is close to the density 
of liquid water. Equation (35) describes the energetics of n(r), and the interface it forms 
is the liquid-vapor interface. 

To the extent that n(r) is a constant equal to the average density of the liquid, Sp{r) 
will have the Gaussian hamiltonian, i.e.. 



where x^(r, r') is the functional inverse of the density-density correlation function, i.e., 
X(r, r') = {dp{r) Sp{r')). The subscript "S" indicates that this hamiltonian applies to 
sma/Z-lcngth scale fluctuations. 

The presence of a large enough solute will force the fluid to be inhomogeneous on 
large length scale. To account for that possibility, the two fields n(r) and Sp{r) must be 
coupled, and the simplest way to do so is with a bi-linear form. Specifically, we take 



Ref. [50]. 



(36) 



p(r) = n(r) + 6p{r) , 



(37) 




(38) 



H[n{r),5p{r)] = Hi^[n{r)] + Hs[Spir)] + Hi[n{r),Sp{r)] 



with 



(39) 
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where -ffnorm['^(r)] ensures that (■^) 

(40) Yl exp{-/3i/[n(r),5p(r)]} = exp{-/3iJL[n(r)]} . 

5p{r) 

The subscript "I" labehng Hi[n(r), dp{r)] stands for interaction, and the symmetric func- 
tion w(r, r') specifies the strength and range of the interaction between the two fields 

While the unperturbed liquid partition sum in Eq.(40) leads back to the Landau- 
Ginzburg description at large-length scales, a non-trivial alteration occurs in the presence 
of imposed inhomogeneity. In the context of possible de-wetting, the most important of 
these alterations comes from excluded volume v due to a solute. This constrains the 
partition sum in a fashion that is compactly described with the functional 



(41) a[p(r)] 



1, when p(r) = for all r G v, 
0, otherwise, 



where v denotes the volume that the solute excludes from the solvent. The volume can be 
complicated, indeed not even contiguous. See for instance the excluded volumes depicted 
in Fig. 9 of Lecture Two. The excluded volume constraint is common to all solutes, 
hydrophobic or hydrophilic. For the latter, a different constraint functional could also 

be employed, one that binds solvent a regions of space adjoining the excluded volumes. 
With the constraints imposed by the solute, the partition sum over small-length scale 
density fluctuations is then 

(42) Yl exp{-/3i?[n(r), 5p(r)]} ^[^(r) + Sp{r)] = exp{-^F,[n(r)]} , 

5p{r) 

where 

(43) Hv[n{v)] = i?L[n(r)] + AH^nir)] 



The alteration to the Landau-Ginzburg hamiltonian, AHy[n{r)], is straightforwardly 
(though tediously) evaluated by carrying out the indicated sum over the Gaussian field 



(^) We leave it as an exercise for the reader to carry out the indicated sum over the Gaussian 
fields to show that i?norm[n(r)] is an irrelevant constant plus 

^ dr y dr' y ds J ds' n(r) M(r, s) x(s, s') u(s', r') n(r') . 
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in Eq. (42) (3). 

In using these formulas to compute numbers, some information about the function 
u(r, r') is required. As it generates a force on the slowly varying field, its functional form 
need not be very specific. It suffices to characterize the function in terms of a mean 
strength and range, i.e., u(r, r') = a(t){\v — r'|), where a is the mean strength with a 
value of the order of k^T, and <;&(r) is normalized with a range of a few A's. The theory 
constructed with Eq.(42) is not terribly sensitive to the specific values of strength and 
range. Physical arguments can be made to estimate their values a priori. Alternatively, 
the parameters can be adjusted so that the theory produces essentially perfect agreement 
between its predictions and the results of a few representative simulation calculations. A 
non-zero value of a is needed to capture the wide breadth of the crossover from small- 
to large-length scale behaviors, e.g., as depicted in Fig. 14. 

A subtlety in these formulas concerns the form of xi'^i r')- When n(r) is not a constant, 
this variance is not simply the liquid phase function given in Eq. (14) of Lecture Two. 
A simple interpolation formula can be used to estimate the variance of the small-length 
scale field, 

(44) x(r,r')«n(r)(5(r-r') + n(r ) [^(Ir - r'|) - 1] n(r'). 

Notice that this formula guarantees that there are no density fluctuations wherever n(r) = 
0. Again, because n(r) is slowly varying, more molecular-scale detail than provided 
reliably by this interpolation formula is unimportant to the evaluation of Affi,[n(r)]. 
This evaluation establishes the following behaviors: 

1. If the excluded volume v is not much larger than the correlation volume of the 
unperturbed liquid, or if it is composed of several distantly separated excluded 
volumes, each one not much larger than a correlation volume, then AiJ„[n(r)] is 
essentially a constant and equal to the solvation free energy of the excluded volume 
as given by Gaussian fluctuation theory, Eq. (12) in Lecture Two. 



(^) We leave it as a second exercise to the reader to show that the result of this calculation is 

A//4n(r)] = -^lndet(x-V2^) 

4-—^ j &r j dr' n(r) -|- J ds J ds' x(r, s) /3u(s, s') n(s') 



Xv (r,r') 



+ 



n(r ) + ds ds x(r i s ) /9w(s , s ) n(s ) 



where Xv ^(r, r') is the inverse of x(i", r') when r and r' are restricted to the volume v. In other 
words, xir^(r, r') satisfies 



dr'x« ^(r, r')x(r', r") = S{r - r"), for r, r" e v. 
Full details on a method of deriving this result can be found in Ref. [51]. 
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2. On the other hand, if v presents a surface of low curvature extending over lengths 
larger than ^, then AH^[n{r)] is no longer constant and is relatively large when 
n(r) ^ for r G To avoid this energetic penalty, probable configurations of the 
slowly varying fields will adjust to make n(r) = within the excluded volume. As 
such, according to Eq. (35), the probability for configurations of the slowly varying 
field will be maximal at configurations with a liquid- vapor-like interface adjacent 
to the excluded volume. 

This latter situation is that of dewetting. Regions of space where dewetting occurs is 
where an excluded volume v causes the average value of n(r) to be zero. It occurs only for 
cases of sufficiently large excluded volumes. The field n(r) governed by Eq. (35) is a con- 
tinuum version of an Ising model or lattice-gas model. We see from the theory sketched 
above that the underlying physics of dewetting is captured by coupling of an Ising-likc 
field to a small-length scale field through excluded volume perturbations. Indeed, with 
a single fixed parameter the strength parameter a, the results of these equations agree 
quantitatively with those of computer simulations, as we have recently shown in detail 
in Ref. [52]. 

3'5. Applications and hydrophobic collapse. - The first general theoretical treatment 
of dewetting in water and its role in hydrophobicity was provided by the Lum-Chandler- 
Weeks (LCW) theory [53]. That theory is the mean-field approximation to what we have 
presented above. In particular, 

(45) (n(r)) « nLCw(r), 
where nLcw(i') is the field that minimizes Hy[n{r)], i.e., 

(46) = 5F4nLcw(r)]/(5nLcw(r). 

The LCW paper [53] presents several illustrative predictions of Eqs. (45) and (46). These 
predictions have motivated many subsequent studies of hydrophobic effects. 

More recent work building from this approach have focused on fluctuations in n(r). 
These fluctuations are important in dynamics. The weight functional is 

(47) P[n(r)] cx exp{-/3F„[n(r)]} , 

which can be sampled by Monte Carlo. The simplest implementation replaces n(r) with 
a binary field on a lattice, prij, where p is the bulk liquid density and rii is either or 
1 depending upon whether ith cell in the lattice is vapor-like or liquid-like, respectively. 
Hi,[n{r)] is then taken to be a lattice-gas hamiltonian, and AHy[n{r)] « Y^^criiVi, where 
c is a positive constant, and Vi is the volume excluded by solutes in cell i. That is. 



(48) 



H[n{r)]K-e'^ niiij + "^{cvi - p fi) rii , 
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where ^ is the chemical potential of the liquid and the primed sum is over nearest neigh- 
bors. The cvi terms account for excluded volume in liquid-like regions, with the free 
energy cost for this excluded to be proportional to the size of that volume. Proportional- 
ity to excluded volume is approximately consistent with small-length scale hydrophobic 
free energies of solvation, as we explained in Lecture Two. 

This version of the theory is particularly easy to implement. With judicious choices of 
lattice spacing and the constant c, most qualitative features of dewctting and hydrophobic 
forces of assembly are captured correctly. The principal feature it fails to capture is the 
slow approach to the macroscopic surface-area scaling illustrated in Fig. 14. Instead of 
a broad crossover from small to large-length scale hydrophobicity, this simplest version 
exhibits a relatively abrupt crossover around 1 nm. This deficiency can be ameliorated, 
as we have detailed in Ref. [52] , but the simplicity of the model described with Eq. (48) 
makes it an attractive choice for easy estimates of roles of hydrophobic forces. 

Figure 22 shows a result obtained from this form of modeling. Depicted are snap shots 
of a trajectory illustrating the collapse of a chain of hard spheres in water. The chain 
follows Newtonian dynamics with friction and random forces reflecting the effects of the 
small-length scale field that has been integrated out. The liquid's dynamics is Monte 
Carlo. Motion of the chain is reflected in changes of w^, and as such the liquid's slowly 
varying fleld couples to the dynamics of the chain. The two move together illustrating 
the collective nature of hydrophobic forces of assembly. 

The intra-chain forces for the chain considered in that flgure are such that the ex- 
tended chain is most stable conflguration in the gas phase. Solvation, in this case hy- 
drophobic forces of assembly, make the globule state the most stable conflgurations in 
the liquid. The half-life of the extended chain in the liquid is about 1 fxs. The pictured 
trajcictory shows parts of a trajectory during the relatively short period of time when the 
chain passes from its metastable extended state to the stable globule state. 

The transition state occurs when a random fluctuation of the chain produces a large 
enough cluster of hard spheres to nucleate a soft liquid-vapor-like interface. At that 
stage, water moves away with relatively low free energy cost and the chain collapses to 
its most stable thermodynamic state. In eff'ect, the reorganization of the chain boils away 
the water. It is a suggestive of the possibility of a nano-scale steam engine, albeit in a 
much more complicated device than a simple hydrophobic chain. It is a possibility that 
we ponder and wonder if it exists in some biological molecular motor. 
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